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CN 1 ABSTRACT 

Aims. The radio-quiet quasar PG 1416-129 (z = 0.129) exhibits atypical optical and X-ray properties. Between 1990 and 2000, in response 
' to its optical continuum decrease, the "classical" broad component of H/? almost completely disappeared, with a factor of 10 decrease in the 
. line flux. In addition, the width of the broad component of the Hy3 line decreased significantly from 4000 km s to 1450 km s . In the X-ray 
O |' band, this object was observed by Ginga in 1988 to have the hardest quasar photon index, with r=l.l±0.1. We present an XMM-Newton/EPIC 
' | observation of PG 1416-129 performed in July 2004. 

>| . Methods. We analyze the time-averaged pn spectrum of this quasar, as well as perform time-resolved spectroscopy. 

' Results. We find that during the present XMM-Newton observation, PG 1416-129 still has a rather hard photon index, both in the soft 
^ (0.2-2keV) and hard (2-12 keV) energy ranges, compared to radio-quiet quasars (BLS1 and NLS1) but compatible with the photon index 
• • ■ value found for radio-loud quasars. This object also shows long-term luminosity variability over 16 years by a factor of three with a variation 
. of photon index from ~1.2 to ~ 1.8. In the soft energy band (0.2-2keV), we found a very weak soft X-ray excess compared to other RQ 
, quasars. The whole time averaged spectrum is fit very well either by X-ray ionized reflection from the accretion disk surface, by a warm 
■ absorber-emitter plus power-law, or by a smeared absorption/emission from a relativistic outflow. While no constant narrow FeK line at 
C3 ' 6.4 keV is observed, we find the possible presence of two non-simultaneous transient iron lines: a redshifted narrow iron line at about 5.5 keV 
(96.4% confidence level according to multi-trial Monte-Carlo simulations) at the beginning of this observation and the appearance of a line at 
6.3-6.4 keV (99. 1 % c.l.) at the end of the observation. These variable lines could be generated by discrete hot-spots on the accretion disk surface. 
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1. Introduction 

PG 1416-129 (z = 0.129) was one of the three bright- 
est low-redshift (z<0 .2) radio-quiet quasars observed during 
the ASCA mission ( Reeves & Turnei 2000l). It has a high 
bolometric luminosity. 1 Operas" 1 ivestereaardl 120021). and 
a black hole mass of about 3-6 x 1O 8 M0 (IVestergaard 



2002; lHao et alJ 120051) . This object has an accretion 



rate with respect to Eddington of about 0.12 (lHao et al 



20051) . It has been associated with the 2MASS object at 
g./7.noo =14 h 19 m 03 s :81, (W)=-13°10'44'/7 dBarkhouse & Hall 
2001). It has strong, broad permitted optical emis sion lines: 
FWHM(H/3)= 61 10km s 1 (IBoroson & Greenlll992l) . but the 
measurement of this width was contaminated by the presence 
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of a very broad component, as shown bv lSulentic et al. (2000), 
who measured FWHM(H/3)= 4000kms _1 . This value of the 
line width led to it being described as a Broad-Line (BL) AGN 
(FWHM(H/3) >2000kms _1 ). However, observations of this 
broad-line quasar have revealed a significant drop-off in its con- 
tinuum luminosi ty by about a factor o f four at 4500 A between 
1990 and 2000 dSulentic et al.ll2000l) . In response to this con- 
tinuum decrease, the "classical" broad component of H/3 has 
almost complet ely disappeared, with a factor of 10 decrease 
in the line flux dSulentic et al. 2000). Additionally, the width 
of the broad component of the H/3 line decreased significantly 
over those 10 years, from 4000km s _1 to 1450kms~'. 
In the X-ray band, PG 1416-129 shows peculiar behav- 
ior compared to other low-redshift quasars. A spectrum ob- 
tained by Ginga in 1988 showed the hardest spectral in- 
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Table 1. Summary of the power-law index and luminosities 
(x 10 44 ergs~') observed by different X-ray satellites in the 
2-10keV energy range over 16 years (5 observations). For 
a direct comparison, the luminosities have been re-calculated 
when needed assuming Ho=75kms _1 Mpc" 1 and qo-0.5. The 
refere nces (c olumn 5) correspond to: RTOO (IReeves & Turner 
2000), LT97 (lLawson & Turned! 19971) . 



Satellites 


Obs date 


r 


Lum 


Ref 




(dd/mm/yyyy) 








Ginga 


02/02/1988 


1.18±0.06 


1.4 


LT97 


Ginga 


21/09/1991 


1.48+0.13 


1.2 


LT97 


ASCA 


29/07/1994 


1.78+0.02 


3.0 


RTOO 


RXTE 


21/08/1998 


1.17+0.15 


1.0 


This work 


XMM-Newton 


14/07/2004 


1.56+0.04 


0.98 


This work 



dex (r= 1.1 + 0.1) of all known qu asars in the 2-20 keV en 



ergy range (Will iams et al.l 119921) . The Compton Gamma- 
ray Observatory/OS SE observations in 1994 showed that 
PG 1416-129 has a very steep spectrum (r=3.2+0.5) above 
50 keV, and consequen tly the most dramatic cut- off between 
50keV and 100 keV dStaubert & Maisackl Il996l). Ho wever , 
from an ASCA observation in 1994. IReeves & Turner (2000) 
found a photon index for this object of F= 1.78+0. 02, which 
is closer to th e values found in mo st radio-quiet low-redshift 
quasars (e.g., iPorquet et al.l l2004al) than the very hard pho- 
ton index measured previously with Ginga. In addition, during 
this ASCA observation, a possible Fe Kor line was detected at 
£=6.54+[J-^keV (E W= 140+75 eV) with an F-test probability 
of 98.3%. This object also shows long-term flux variability as 
reported in Table [TJ In this paper we report the analysis of the 
XMM-Newton observation of PG 1416-129 performed in July 
2004. In addition to the previous X-ray observations already 
published and the present XMM-Newton observation, we also 
show the currently unpublished 1998 RXTE data. Details of the 
RXTE data analysis are reported in AppendixlAl 



2. XMM-Newton observation and data reduction 

XMM-Newton observed PG 1416-129 on July 14, 2004 
(ID 0203770201; orb it 842; exposure time ~50ks). The 
EPIC-MOS cameras dTurner et al.l200lh operate d in the Small 



Wind ow mode, while the EPIC-pn camera ( S tra der et al 



2001) was operated in the standard Full Frame mode. The data 
were re-processed and cleaned (net pn time exposure ~45ks) 
using the XMM-Newton SAS version 6.5.0 (Science Analysis 
Software) package. Since the effect of pile-up was negligible, 
X-ray events corresponding to patterns 0-12 and 0-4 events 
(single and double pixels) were selected for the MOS and pn, 
respectively. Only good X-ray events (with FLAG=0) were in- 
cluded. The low-energy cutoff was set to 300 and 200 eV for 
MOS and pn, respectively. In this article we present the pn 
results only, since the pn-CCD have a better sensitivity over 
a broader energy range (0.2-12keV) compared to the MOS 
CCDs. However, we verified that similar results were obtained 
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Fig. 1. The pn spectrum of PG 1416-129 (in the observed 
frame). A power-law has been fit to the 2-5 keV data and ex- 
trapolated to lower and higher energies. A weak soft X-ray pos- 
itive residual is seen between 0.3-1 keV, as well as the presence 
of a Fe Ka complex between about 4.8-5.7 keV (~5. 4-6.4 keV 
in the quasar frame). For presentation only, the data have been 
re-binned into groups of 10 bins, after a grouping of a mini- 
mum of 20 counts per bin is used for the fit. 

using the MOS data. The pn source spectra were extracted us- 
ing a circular region of diameter 30" centered on the source 
position. PG 1416-129 is by far the brightest X-ray source in 
this 30' EPIC field-of-view. Background spectra were taken 
from a box region (7.3' x 4') on the same CCD as the source 
(excludin g X-ray point sources). The xspec vll.3 software 
package dArnaud 1996) was used for spectral analysis of the 
background-subtracted spectrum, using the response matrices 
and ancillary files derived from the SAS tasks rmfgen and ar- 
fgen. The pn spectrum was binned to give a minimum of 50 
counts per bin. 

The signal to noise ratio was not sufficient for reliable RGS 
data analysis, and the OM was blocked during this observation. 

3. XMM-Newton spectral analysis 

In all fits, we included both the Galactic column density 
(/V^ al =6.80 xl0 20 cirT 2 , obtained from the colde ns program 
using the compilations of lDickev & Lockmanll 990). and a pos- 
sible additional absorption component located at the quasar 
redshift (N 1 ^). Note that all fit parameters are given in 
the quasar's rest frame, with values of Ho=75kms~ 1 Mpc -1 
and qo=0.5 assumed throughout. The errors quoted corre- 
spond to 90% confidence ranges for one interesting parame- 
ter (A y 2 -2.71). Abundances are those of lAnders & Grevesse 
dl989l) . In the following, we use the updated cross-sections for 
X-ra y absorption by the interstellar m edium (tbabs in XSPEC) 
from lWilms. Allen. & McCravl d2000t) . 



3.1. The time averaged spectrum 

We found an average pn count rate in the 0.2-12 keV energy 
range of 1.36+0.01 cts s _1 . First, fitting the overall 0.2-12keV 
pn spectrum with a single absorbed power-law model, a mod- 
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erately acceptable fit was found Cfcd d = 1.18 for 979 d.o.f.), 
with r= 1.75 ±0.01. To characterize the hard X-ray continuum, 
we fit an absorbed power-law model over the 2-5 keV energy 
range, where the spectrum should be relatively unaffected by 
the presence of a broad soft excess, a Warm Absorber-Emitter 
medium, an FeKa emission line, and a contribution above 
8keV from a high energy Compton reflection hump. In this 
energy range, the data are well fit by a single power-law model 
with T= 1.61 ±0.06 Ct 2 /d.o.f.=346.2/377). This power-law in- 
dex is similar to values found over the same energy range 
for radio -loud quasars (<r= 1.74> with a standard deviation 
of 0.03: IPorquet et al" 2004aj), but harder than the values com- 
monly seen in radio-quiet quasars: <F=1.90> (with a standard 
deviation of 0.27) and <F=2.37> (with a standard deviation 
of 0.11) fo r Broad Line quasars and Narrow Line quasars, re- 
spectively dPorquet et alj2004al) . A similar conclusion is found 
when comparing with the average power law index measured 
in radio-quiet qu asars above 2 keV [< r2-iokeV>= 2 with a dis- 
persion of 0.25. iGeoree et alJ l2000t <r 2 _i 2ke v>=1.89±0.11, 
Piconcelli et alj|2005l) . 

Figure Q] displays the spectrum when the power-law model 
is extrapolated over the 0.2-12 keV broad band energy range. A 
weak positive residual is seen below 1 keV due to the presence 
of a soft X-ray excess. In addition, a Fe Ka complex between 
about 4.8-5.7 keV (~5 .4-6.4 keV in the quasar frame) seems to 
be present. 



3.1 .1 . The soft spectrum: 0.2-2 keV 

First we fit the 0.2-2 keV energy range with a single 
absorbed power-law, and found a good fit with NS l - 



2.1±0.7 x 10 zu cnr 2 and a photon index of F=1.96±0.05 
(X 2 /d.o.f. =364.9/355). This value is compatible within the er- 
ror bars with T=2.2±0.2 i nferred from a ROSAT P SPC obser- 
vation performed in 1992 (Ide Kool & Meurslll994l) . This value 
is marginally smaller than the one found in general in radio- 
quiet quasars (Por quet et al.ll2004at <ro.3_2kev>=2.43 with a 
dispersion of 0.35). The inferred 0.2-2 keV unabsorbed flux 



is 2.9x lO-^ergcnrV 1 (L=l.l x lO^ergs- 1 ). However, 



adding a narrow (cr=10eV) Gaussian line, we found a sig- 
nificant improvement to the fit, with ^ 2 /d.o.f. =339.4/353 (as 
a comparison, Monte Carlo simulations, as described in 93.21 
give >99.9%). The energy of the line is 556±15eV and its 
equivalent width (EW) is 15.9^ eV. These parameters are 
consistent with the position of the He-like O vn triplet between 
561-574 eV. We also fit the data assuming two absorption 
edges, and we again found a significant improvement to the 
fit, with ^ 2 /d.o.f =339.9/351 (F-test >99.99%). The energies 
(and optical depths, t) of the two edges are E-Q.llA^^keW 
(r=0.16±0.06) and £=0.953!°;^ keV (t=0.10±0.05) respec- 
tively. The energy of the first edge is consistent with the 
rest-frame energy of the Ovu edge (0.740 keV). The second 
edge is not compatible with the rest-frame of the O vm edge 
(0.871 keV) or of Ne ix (1.196keV), but could be one of these 
ions with a velocity shift or perhaps a blend of the two edges. 
It also could be due the presence of absorption lines, for ex- 
ample, the Fe L complex seen in absorption in several AGN 



(e.g. JSako etaDl200lh . Unfortunately the present RGS data do 
not have good enough S/N to determine the physical nature 
of this feature. However, as we will discuss below, the whole 
average spectrum could be equally well explained either by X- 
ray ionized reflection on the accretion disk surface, by a warm 
absorber-emitter plus power-law model, or by a smeared ab- 
sorption/emission from a high velocity outflow. 

3.1 .2. The Fe Ka complex 

As shown in Fig.Q] there is the presence of a FeKa: complex 
between 5.4 and 6.4 keV in the quasar frame. Therefore, we fit 
the 2-12 keV energy range with an absorbed power-law con- 
tinuum and three possible profiles for the iron line; either (i) 
a simple Gaussian (zgauss in xspec), (ii) a line profile corre- 
sponding to emission from a relativis tic accretion dis k around 
a non-rotating black hole (diskline; iFabian et al . 1989), or (iii) 
emission from around a maximally rotating black hole (laor; 
Laoi 199 lh . The spectral fits were improved by the addition 
of the iron line, when modeled either by a Gaussian line or 
by a relativistic accretion disk-line profile, with a significance 
greater than or equal to 99.3% from an F-test (see Table |5J. 
The line, if modeled by a Gaussian, appears to be more likely 
narrow (cr <0. 1 1 ke V) and produced by neutral to moderately 
ionized iron, with £'=6.34±0.05keV. If the line is due to a 
relativistic accretion disk, the inclination of the disk is con- 
strained to be less than about 20°. The line could be associated 
with X-ray fluorescence from the molecular torus. However, 
we will see below with time-resolved spectroscopy that the sit- 
uation is more complex, with the possible presence of two non- 
simultaneous transient iron lines at ~5.5ks and ~6. 3-6.4 keV 
appearing during the first 5.5 ks and last 16 ks of the observa- 
tion, respectively. 



3.1.3. The broad-band spectrum: 0.2-1 2 keV 

The 0.2-12 keV energy range is well-fit by the com- 
bination of a blackbody and a power-law component. 
We find A^ nt < 2.4+J^xlO 20 cirT 2 , /tr bb =149! 2 ° eV and 
r=1.62+^f| Cr 2 /d.o.f =1008.7/977). A broken power-law 
model also gave a good fit, with N^= 2.9+jjXlO 20 cirT 2 , 
F soft =2.04±»»*, F hard = 1.59^|» and £ break =L64±°- 2 !; keV 
Ot 2 /d.o.f =995.8/977). 

In order to obtain a more physical representation of the 
soft excess, we have also tested multi-temperature disc mod- 
els, which may be expected if the soft X-ray excess origi- 
nates via thermal emission from the inner accretion disc in 
PG 1416-129. The diskbb (non-relativistic) and diskpn (rel- 
ativistic) models were used, together with a power-law to 
model the hard X-ray emission above 2 keV. Equally good fits 
are obtained for both models. For the first model, we found 
£r diskbb =211+^eV and F=1.60±0.03 Ct 2 /d.o.f. =996.0/977), 
and similar parameters were found for the second model: 
£7d iskpn = 200+geV and F=1.60±0.03 Ct 2 /d.o.f. =996.3/977). 
The inner disc temperatures obtained through either of 
these models are much larger than the temperature at 3Rs 
(/?s=2GM/c 2 ) expected from a standard steady state a thin ac- 
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Table 2. Best-fit spectral parameters of the time-averaged pn spectrum in the 2-12 keV energy range with an absorbed (Galactic, 
A^H=6.8xl0 20 cm 2 ) power-law (PL) component plus a line profile: zg auss: Gaussian prof ile; and diskline and laor: line profile 
emitted by a relativistic accretion disk for a non-rotating black hole dFabian et al. 1989 ) and a maximally rotating black hole 
dLaorll 19911) . respectively. The line fluxes are expressed in 10 6 photons cm 2 s 1 . We assume an emissivity index q equal to -2. 
(a): R m =6R g and/? out =1000fi g . (b): R in =l.235R g and R ont =400R g . 



Model 


r 


E (keV) 


Line parameters 
cr (keV) 6 (deg) Flux 


EW (eV) 


x 2 /d.o.f. 


F-test 


PL 


1.56±0.04 












634.0/621 




PL + zgauss 


1.57±0.04 


6.4 (f) 


0.01 (f) 




2.2±1.2 


54±29 


624.5/620 


99.78% 


PL + zgauss 


1.57±0.04 


6.33±0.05 


<0.10 




2.9+1 .4 


69+36 


619.6/618 


99.74% 


PL + diskline'" 1 


1.58±0.04 


6.40±0.05 




<22 


3.9+1.7 


107±46 


619.3/618 


99.77% 


PL + LAOR (6) 


1.58±0.04 


6 40 +0 05 




<22 


4.4±2.1 


119** 


621.7/618 


99.30% 



cretion disc, i.e. 8-10 eV, dShakura & Sunvaev 1 19731) . assum- 
ing a black hole mass of about 3-6xl0 8 Mq. This result is con- 
sistent with the temperature found in other low-redshift AGN, 
with a 1 00-200 eV temperature component being too hot to 
be the direct thermal disk emission as inferred from the soft 



X-ray spectra ( 


Gierlinski & DondEooH IPorquet et al.1 l2004ji 


Piconcelli et al. 


2005). 



The whole spectrum could be explained with a warm 
absorber/emitter combined with a power law continuum. 
Therefore, we used a grid of mod els generated by the xstar 



(version 211) photoionization code (IKallman et al.ll 19961) to fit 



the EPIC-pn data. A one zone model was adopted (covering 
the range of parameters l<log£ <5 and 2 1 <log /V™ arm <24). 
Solar elemental abundances were assumed. A low turbulence 
velocity of 200km s _1 was used, and the covering factor was 
fixed to unity. The warm absorber column density (/V™ 1 ™) was 
fixed to the lowest value in the xstar grid, 10 21 cirT 2 , while 
the ionization parameter (£) was left free. A good fit was 
obtained Ct 2 /d.o.f.=1077.6/975) with JVjf <1.6xl0 19 cirT 2 , 
r=1.69+0.02, log £=1.90+°^ ergons" 1 . In this model, the 
photoionized emission from the warm absorber represents 
about 7% of the total flux in the 0.2-2 keV energy range. 

We then tested the possibility that the whole spectrum 
could be explained by X-ray ionized reflectio n on the accre- 
tion d i sk surface (reflion m odel in XSPEC; Ross & Fabianl 
2005). Ross & Fabianl (2005) have shown that this model can 
explain the overall spectrum from 0.2-12 keV in several AGN, 
with out invoking ad-hoc components such as a blackbody (see 
also lCrummv et alj|2005l) . One signature of ionized reflection 
is the soft excess emission which occurs in the 0.2-2 keV band 
due to lines and bremsstrahlung emission from the hot surface 
layers. It creates a bump in the relativistically blurred spec- 
trum which, when convolved with the XMM-Newton pn re- 
sponse matrix, is well fit by a blackbod y of temperature 150 eV. 
The i ron abundance (A Fe ) is from iMorrison & McCammon 
( 11983b in this model. The reflion model covers a large range 
of parameters: 30 < £ < 10000 ergcms -1 (ionization pa- 
rameter), 0.1 < Ap e < 10.0 (relative iron abundance), and 
1.0 < r < 3.0. We also took into account the relativistic ef- 



fect due to relativistic motion in the inner part of the ac- 
cretion disc, by blurring the spectrum with a Laor line pro- 
file (kdblur; IRoss & Fabianl 120051: ICrummv et ail 120051) . In 
the blurred reflection model, the physical lower limit of the 
inner radius (Rm) was set at 1.235 R s , i.e. corresponding to 
the innermost stable orbit around a maximally rotating black 
hole. We fixed the disc emissivity index to -3. We obtained 
a very good fit 0f 2 /d.o.f.=980.9/974) to the data with the fol- 
lowing parameters: /Vjf=3.6 + 1.3 x 10 20 cirT 2 , R in ^3.l^R g 
(R g =GM/c 2 ), <9=21.1+|°- 2 °, A Fe =0.92+°^ r=1.78+0.04,'and 
^=38.0^ 2 ' 8 ergcm _1 s -1 (see Fig. |2). No additional absorp- 
tion or emission features in the soft band were required. As 
shown in Fig. [2] the reflection emission is significant in the 




0.2 



0.5 



1 2 
Energy (keV) 



10 



Fig. 2. The time-averaged pn spectrum of PG 1416-129 (in the 
observed frame ). A relativistically bl urred photo-ionized disc 
model (reflion: IRoss & Fabianl 12 0051) has been fit, and which 
provides a good representation of the whole energy band. A pri- 
mary power-law component is taken into account in the model. 
Both component models are displayed without blurring in the 
upper panel. See text for the values of the fit parameters. For 
presentation only, the data have been re-binned into groups of 
10 bins, after group of a minimum of 20 counts per bin is used 
for the fit. 
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Table 3. Best-fit spectral parameters for the three pn sub- 
spectra in the 2-10keV energy range. Column (1): number of 
the corresponding sub-spectrum (the extraction time range is 
reported on Fig. [3]). Column (2): model used. PL: power-law; 
zgauss: Gaussian line (<x fixed to 10 eV); AN: narrow diskline 
annulus (AR=lR g ; other diskline parameters are discussed in 
the text, see 33.21 ). Column (3): power-law index. Column (4): 



energy (/keT 
8 



unabsorbed 2-10keV flux (xlO 12 erg cm 1 s '). Column (5): 
line parameters. 



n° Model 


r 


runabs 
r 2-l0 


Em* (keV) 
EW (eV) 


*7a.o.f. 


1 PL 

PL+zgauss 


1.61±0.11 
1.64±0.11 


3.2+0.5 


5.47+0.04 
194+89 


83.8/97 
71.0/95 


2 PL 


1.51+0.05 


3.3+0.2 




98.8/97 


3 PL 

PL+zgauss 

PL+AN 


1.61±0.06 
1.62+0.06 

1.64+0.06 


3.1+0.2 


6.35+0.05 
128+53 
6.40 (f) 
287+100 


117.9/97 
102.5/95 

95.0/93 



very soft part of the spectrum, especially near 0.4-0.5 keV (ob- 
served frame), where a feature at 556 e V has been detected us- 
ing a simple power-law plus a narrow Gaussian (see §13.1.11 1. 
The fraction of the reflection emission to the total emission is 
only about 21% in the 0.2-12keV (similar values are found in 
the 0.2-2 keV and 2-10keV energy ranges). The inner radius 
found here is compatible with a rotating black hole. The fit pa- 
rameters found here are compatibl e to those determined in the 
large sample of type I AGN from ICrummv et all d2005l) . The 
power-law index is in the low range value of the sample but not 
as low than that found for Mrk 359 (r= 1 .49+0.06), which is an 
extreme NLS 1 . 

An alternative possibility is that the soft X-ray excess is 
an artifact of soft X-ray absorption and emission produced 
by a h igh velocity outflow, as show n by iGierlihski & Done 
d2004 and ISchurch & Pond d2006l) . who model the X-ray 
spectra of other PG quasars observed by XMM-Newton. 
Therefore we fit the data with smeared ionized absorption 
and emission, using the same grid of models generated by the 
xstar photoionization code used for the warm absorber-emitter 
model used above. We fixed N™ m to 10 21 cm~ 2 . We obtained 
a very good fit to the data Of 2 /d.o.f. =990.8/975), with a 
smearing velocity of 0.18+0.02c, A^?=9.0 + Hxl0 19 
log £=2.31+° ^ erg cm" 1 s ' 



and T=1.59!° : ° 2 . 



-3.<r 



cm 



As we will see in the next section, an additional physical 
process, such as hot spots due to X-ray flares, is required to 
explain the two possible transient Fe Ka lines observed ( jH3.2l i. 




20000 30000 40000 

time (/s) 




Fig. 3. The signal-to-noise deviations (in cr) of the Fe Ka line 
above the power-law continuum in the energy-time plane. The 
energy axis is given in the observed frame. The colour-scale 
represents excess signal/noise in the line counts above the fitted 
continuum. Energy and time are oversampled by a factor of 
10. Data are top-hat smoothed by 10 ks in time and Gaussian 
smoothed with FWHM =0.14keV in energy. The regions 1- 
3 correspond to the time selections used for the time-resolved 
spectroscopic analysis. See text for explanation. 

3.2. XMM-Newton time-resolved spectroscopy 

To test possible rapid energy shifts of the iron line, we have 
created X-ray intensity maps in the energy-time plane using 
the pn data. Photons from the source cell were accumulated 
in pixels in the energy-time plane. The pixel distribution was 
smoothed in energy by the instrumental resolution, using a 
Gaussian of 140 eV (appropriate for the single plus double 
events with the latest calibration), and smoothed in time us- 
ing a top-hat function of width 10 ks. Each time-slice was 
background-corrected by subtracting a time-dependent back- 
ground spectrum measured in an off-source region on the same 
detector chip as the source. The source continuum was mod- 
eled as an absorbed power-law, with variable amplitude and 
slope but time-invariant absorption column density. This con- 
tinuum was subtracted, leaving positive and negative residuals 
that comprise noise plus any emission or absorption compo- 
nents on top of the cont inuum. More info rmation about this 
method can be found in Tur ner et al ] d2006h . The "signal-to- 
noise" (S/N) map presented in Fig|3]is the ratio of the fluctua- 
tion amplitude to the calculated noise. 

The highest S/N observed is about 4.5. We observe the pres- 
ence of an excess with S/N >3.5 at ~5.5keV in the quasar 
frame at the beginning of this observation. Then this feature 
disappears. At the end of this observation, there is the pres- 
ence of an excess near 6.3 keV (in the quasar frame) and then 
at lower significance (S/N <3.5) a hint of a shift to 6.4keV of 
the centroid of the line. In order to characterize these features, 
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Fig. 4. The three pn time-resolved spectroscopy spectra (in the 
observed frame) obtained by splitting the average spectra ac- 
cording to the time regions shown in Fig. 3. The three sub- 
spectra have been fit with a power-law continuum model. Here, 
5.67 keV in the observed frame corresponds to 6.4keV in the 
rest frame. Notice the residuals present at 4.8 keV and 5.6 keV 
in segments 1 (upper panel) and 3 (lower panel), respectively. 



which are most probably associated with Fe Ka line emission, 
we have split the observation into three sub-spectra (the corre- 
sponding time ranges are illustrated on Fig.O: 
- 1 . This corresponds to the presence of the 5.5 keV feature with 
S/N > 3.5 (from t=0 ks to t=5.5 ks); 

- 2. This corresponds to the absence of any significant long- 
lived excess feature (from t=5.5ks to 32 ks); 

- 3. This corresponds to the appearance and presence of the 
feature near 6.3-6.4keV. 

First, we rebinned the sub-spectra, grouping every 80 eV 
(16 channels). This was chosen to sample the resolution of 
the pn detector (i.e. about 2 spectral bins per FWHM res- 
olution element). We then fit them with a power-law con- 
tinuum over the 2-10keV energy range. The parameter fits 
are reported in Table [3] and the spectral fits are displayed in 
Figure |4] For sub-spectra 1 and 3, an emission feature is ob- 
served. For sub-spectrum 1, the addition of a narrow Gaussian 
line (width fixed to 10 eV) with a centroid energy of ~5.5 keV 
improved the fit. Allowing the line width to vary, we found 
an upper limit of 94 eV. For sub-spectrum 3, the addition of 
a narrow Gaussian (width fixed to lOeV) with a centroid en- 
ergy of about 6.35 keV also improved the fit. If we allow the 
line width to be a free parameter, the upper limit is rather 
large (<x <0.4keV) due to the presence of the red wing-like 
feature seen below 6.4 keV (see Fig. HJ). A narrow disc annu- 
lus at 6.4keV (AR-1 R g , emissivity index fixed to -2) pro- 
vides a superior fit. The best fit parameters are T= 1 .64+0.06, 
Rin=34tlR s , 6>=17+3° and £W=287+100eV (Fig.0. Letting 
the value of the accretion disc outer radius free to vary, we find 
an upper limit of 84 R g , corresponding to a AT? of the annulus 
value between 1 and 58 R g (taking into account the lower limit 
on R in ). 
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Fig. 5. The pn sub-spectrum 3 of PG 1416-129 (in the observed 
frame). A power law continuum plus a diskline model have 
been fit. A narrow annulus is assumed with AR=1 R s . See text 
for the best fit parameter values. 



As 



Porquet et al. 



discussed by 
2004bl) . 



Protassov et al 



(|2002h (see also 
the F-test can overestimate the true 
detection significance of a line if there was no previous 
expectation of the line. One must account for the number 
of resolution elements N over the bandpass of interest. The 
probability of detecting a feature at any energy in this range is 
obtained using (P\) N , where P\ is the probability of finding 
a feature at an expected line energy (equal to one minus the 
null hypothesis probability). Assuming 21 resolution elements 
between 4-7 keV, the significance probabilities for the 5.5 and 
6.3 keV lines become 98.1% and 99.2%, respectively. We also 
carried out a more rigorous test of the significance of the lines 
using Monte Carlo simulations to assess the likel ihood that the 
lines c ould be an artefact of photon noise (e.g., IPorquet et al 



2004b). For our null hypothesis, we assumed that the spectrum 
is simply an absorbed power-law continuum, with the same 
parameters as the absorbed power-law model fit to the real 
data. We used the xspec fakeit none command to create a fake 
spectrum, and then, to account for the uncertainty in the null 
hypothesis model itself, we refit the null hypothesis model and 
ran fakeit none a second time on the "re-fit" null hypothesis 
model. Both fakeit none commands were run with the photon 
statistics appropriate for a 5.5 ksec exposure (segment 1) or 16 
ksec exposure (segment 3). We grouped each faked spectrum, 
rebinning every 80 eV. We then added a narrow (<x = 10 eV) 
Gaussian to the fit for the 5.5 keV line, or a diskline annulus, 
with the values of the inner and outer radii, emissivity index, 
and inclination angle fixed to that of the observed spectrum, 
for the 6.3 keV line. We searched over 4.0 to 7.0 keV (observed 
frame) in increments of 0.1 keV, fitting separately each time 
to ensure that the lowest value of x 2 was found. We compared 
the minimum value of y 2 obtained with the corresponding 
X 2 of the null hypothesis fit. This process was repeated 1000 
times, yielding a cumulative frequency distribution of the 
simulated A% 2 expected for a blind line search. The inferred 
probabilities to reject the null hypothesis are 98.8% and 99.7% 
for the 5.5 keV and 6.3 keV lines, respectively. However, since 
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we have divided the whole observation into three sub-spectra, 
the null hypothesis probabilities have to be multiplied by a 
factor of three, leading to 96.4% and 99.1% for the detection 
confidence level of these transient lines. While the detection of 
the 5.4 keV feature is statistically marginal, it should be noted 
that the 6.3 keV feature is observed at an energy where a iron 
Ka diskline would be expected a priori. 

The three sub-spectra display no significant flux variations: 
is consistent with being constant. As far as the con- 
tinuum photon index is concerned, sub-spectra 1 and 3 have 
consistent values, but sub-spectrum 2 is slightly flatter; the 
photon index when the diskline annulus model is applied to 
sub-spectrum 3 is not consistent at 90% confidence with that 
obtained in sub-spectrum 2, where no emission line was found. 
However, we can use Monte Carlo simulations to quantify 
the expected spread in the measured photon index, given an 
assumed intrinsic constant photon index. Specifically, the 
null hypothesis model in the Monte Carlo procedure above 
assumed a single value for the intrinsic photon index; using 
those same simulated spectra, we find a standard deviation 
in the measured photon index of 0.165 (due to the moderate 
signal/noise). This is higher than the maximum AT of 0.13 
between the best-fit photon indices for sub-spectra 3 and 2. 
Therefore, there is no evidence for the appearance of the line 
(or lack thereof) to be correlated with values of the photon 
index. 



4. Summary and discussion 

In this section, we summarize the main results of our anal- 
ysis of the July 2004 XMM-Newton/EPIC observation of 
PG 1416-129. 

• A X-ray flat RQQ object. 

With the present 2004 XMM-Newton observation, we found 
that its 2-5 keV power-law index (F ~1.6) is similar to values 
found f or radio-loud objects ( T- 1 .74 with a standard deviation 
of 0.03; |Porquet et alj |2004a). but by far much harder than the 
value commonly seen in radio-quiet quasars: T=1.90 (with a 
standard deviation of 0.27) and F=2.37 (with a standard devia- 
tion of 0.1 1) f or Broad Line quasa rs and Narrow Line quasars, 
respectively ( Porquet et al. 2004ah . This result confirms the 
behavior seen in previous observations, with the exception of 
the 1994 ASCA observation, as reported in Table Q] 



• IsPG1416-129aNLSl ? 

PG 1416-129 could be associated with the NLS1 class of 
objects, given the va lue of FWHM(H/3) observed in 2000 by 
Sulenficetai]d2000h . 

NLSls are defined according to the following characteristics: 
- FWHM(Hp) less than 2000 km s~'; strong optical Fen 
multip l ets; \Q Iif | /15007 A/H6 n < 3 dOsterbrock & Pogge 



19851; iGoodrichl Il989h . The 2000 optical spectrum of 



PG1416-129, however, show ed only a weak Fen multiplet, 
with [O in] A5007A/U/3 -13.9 dSulenfic et alfeOOOh. 

' 19921: 



steep, soft X -ray excess dPuchnarewicz 



et alJ 

Boile r et al. Il996l) : rapid soft/hard X-ray variability dLeighlvl 
19991) . During the 2004 XMM-Newton observation, however, 



neither significant X-ray variability nor a steep soft excess 
were observed. 

- less massive black holes with higher Eddington ratios, 
suggestin g that they might be in the early stage of A GN 
evolution dGrupdl 19961; lMafhurj|2000l; iBian & Zhaoll2003l) . As 
reported in the introduction, however, this object has a black 
hole mass of about 3-6xl0 8 M Q , which is higher than the 
mass found in typical NLSls. Furthermore, its accretion rate 
with respect to Eddington is only 0. 12. 

Therefore it appears that PG 1416-129 is not a genuine NLS1 
and instead the narrow Hf3 width is caused by the variations in 
the broad component of this line. 

• The soft X-ray spectrum. 

In addition to its flat 2-10keV photon index, this object 
shows a much weaker soft excess than generally observed 
for Broad-Line quasars and Narrow-Line quasars. Some 
additional features seem to be present in the 0.2-2 keV energy 
range. For example, the addition of a Gaussian emission line at 
556+15 eV improved the fit very significantly, with an F-test 
significance of >99.99%. However, we have obtained a fit with 
the same statistical quality by adding two absorption edges. 
Unfortunately the present RGS data do not have sufficient 
signal to noise ratio in order to determine the physical nature 
of this feature, if real. The whole average spectrum is very well 
fit by a model with X-ray ionized reflection on the accretion 
disk surface, with no additional soft X-ray absorption and/or 
emission features required. However other models, such as 
absorption and/or re-emission from a high velocity outflow 
(or simple photoionized emission from distant matter) could 
also explain the weak soft ex cess below 1 keV, as discus sed in 
detail for other PG quasars bv lGierlihski & Pond (12004b . 



• Discovery of two possible transient Fe Ka lines 
In the time-averaged spectrum, a FeKc complex is seen be- 
tween 5.5 and 6.4keV, with significant residuals at about 
6.4 keV. If the line is due to a relativistic accretion disk, the 
inclination of the disk is constrained to be less than about 20°. 
If the line is variable, this would rule out X-ray fluorescence 
from distant matter such as the molecular torus. Indeed, time- 
resolved spectroscopy has revealed that the emission complex 
is more likely due to the presence of two transient iron lines. 
The first iron line has a redshifted energy centered at ~5.5 keV 
and was present during the first ~5.5ks (EW ~167eV) of 
this observation with a 96.4% confidence level according to 
Monte-Carlo simulations (see 33.21 ); after 5.5 ksec, it disap- 
pears sharply. The second, more significant transient iron line 
appears about 32 ks after the beginning of this observation with 
a centroid energy at ~6.3keV (EW ~250eV) and a duration 
of at least ~16ks with a 99.1% confidence level according to 
Monte-Carlo simulations. There is a hint of a slight shift in its 
centroid energy to ~6.4keV during the last lOks of the obser- 
vation. Other narrow, highly-redshifted iron lines (~4.5-6 keV) 
in AGN have already been re ported: e.g., NGC35 16 (Turner 
et al. 2002) ESO 113-G010 JPorauet et al.1 l2004bh. J 
dTurner et alJl2004, AX J0447-0627 dDella Ceca etal 



2005) 



and PG 1425+267 dMiniutti & Fabianl 120061) . Localized spots 
or narrow annuli which occur on the surface of an accretion 
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disk following its illumination by flares have b e en proposed 
to explain these features ( e .g., Navakshinl 2000t Turner et al 

— 3 



20021 iDovciak et alj |2004 Ipechacek et al.1 l2005h . Such lines 



may arise in the accretion disk, sporadically illuminated by 
"hotspots," as regions of magnetic reconnection illuminate very 
small areas of the accretion flow. Presumably the hotspot co- 
rotates with the disk. While the small line width implies that 
the emitting region must be small and detected only during a 
fraction of the whole orbit, large spots, meanwhile, would pro- 
duce broader features and they would be prone to rapid de- 
struction. Here the features appear to be narrow, however their 
widths are not strongly constrained. Indeed, an annulus of up to 
587? g in size is allowed in segment 3 of the observation. If the 
transient lines are due to hotspots co-rotating with the disk, a 
lower limit of 16 ks on the lifetime of the transient 6.3 keV iron 
line feature would correspond to at least 80% and 2% of the 
whole orbit at R=l .235 R g and R-2QR g , respectively (assum- 
ing a black hole mass of 3 xlO 8 M0). The large EW of both 
transient features are consistent with the X-ray emission being 
located in a relatively compact region, e.g., from X-ray flares. 
More complex geometries of the X-ray emitting region can be 
considere d, e.g., spiral waves propagating across the accr etion 



disc (e.g.. iKaras et al.l 12001: Fuk umura & Tsurutal 12004). and 



density inhomogeneities, which c an indicate turbulence dr iven 
by magnetorotational instability (Balbus & Hawlevlll99ll) . or 
ph oton bubble inst abilities dAronsI 1 1992^ 7 As demonstrated 
by iBallantvne et al] d2004l 120051) . Photon bubble instabilities 
(which result on time-varying density homogeneities on scales 
smaller than the disk thickness) can explain rapid variability of 
emission lines, with EW variations as large as ~100 eV that are 
independent of the illuminating continuum, as observed dur- 
ing the present observation. In addition, scenarios involving 
non-flat ac cretion disks, such as warped/corrugated accretion 



disks (e.g.. lFabianet al. 2002; Miller & Homan 2005), can be 



invoked to explain such features. However, as pointed out by 
Fabian et al.l d2002|) , even a small corrugation of the accretion 
disc may result in T >2 and a strong reflection component in 
the observed spectrum, as observed in several NLS1 objects 
(e.g, 1H 0707-495). In the case of the present observation of 
PG 1416-129 such a scenario appears unlikely since neither 
a steep spectrum (F ~1.6) nor a strong reflection component 
(the fraction of the 0.2-12 keV reflection emission is only about 
21% of the total emission, see 83. 1.3b are observed. 
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Appendix A: RXTE data 

PG 1416-129 was observed by RXTE on 1998 Aug 21 from 
06:25-10:42 UT, 1998 Aug 22 from 06:54-08:06 UT, and 
1998 Aug 23 from 01:37-07:32 UT, for a total of 41.1 ksec. 



Data were taken u sing RXTE's proportional counter array 
(PCA;|Swankjj998), which consists of five identical collimated 
proportional counter units (PCUs). Data were collected from 
PCUs 0, 1, and 2, and reduced using standard extraction meth- 
ods and FTOOLS v5.3.1 software. Data were rejected if they 
were gathered less than 10° from the Earth's limb, if they were 
obtained within 30 min after the satellite's passage through the 
South Atlantic Anomaly, if ELECTRON0 > 0. 1 , or if the satel- 
lite's pointing offset was greater than (f.02. 
As the PCA has no simultaneous background monitoring ca- 
pability, background data were estimated by using pcaback- 
est v3.0 to generate model files based on the particle-induced 
background, SAA activity, and the diffuse X-ray background. 
The 'L7-240' background models appropriate for faint sources 
were used. This background subtraction is the dominant source 
of systematic error in RXTE AGN monitoring data (e.g., 
Edelson & Nandralll999h . Counts were extracted only from the 
topmost PCU layer to maximize the signal-to-noise ratio. This 
reduction yielded a total of 22.0 ksec of good exposure time. 
Response matrices were generated using pcarsp v. 8.0. For the 
purposes of spectral analysis, data below 3.0 keV were dis- 
carded in order to disregard PCA calibration uncertainties be- 
low this energy. The source was not well detected by the PCA 
above ~14 keV; photons above this energy were discarded. 

Fitting these RXTE data with an absorbed (A^ al ) power- 
law continuum, we found a good fit with r=1.17+0.15 
0f 2 /d.o.f.=7.6/25). The inferred 2-10keV luminosity is 1.0 x 
ergs~ 
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